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Through the last years, several types of numerical and combinatorial optimization algorithms 
have been used as useful tools to minimize functional forms. Generally, when those forms are 
non-linear or occur in problems without a specific optimization method, stochastic methods based 
on search algorithms have shown good results due to its smaller susceptibility to be trapped in a 
local minimum. Besides that, they can easily be implemented to work with other techniques, in 
this class of algorithms, the genetic ones have received special attention because they are a robust 
optimization tool. An algorithm can be named genetic when it uses some kind of codification to 
transform a set of possible solutions of a given problem in a population that will evolve subject to 
operators inspired, or not, by mechanisms of natural selection. In other words, they work with a 
population of solutions to obtain better solutions in the next generation. To do this, they use only 
information of cost and prize. In this work, we propose a genetic algorithm optimization technique 
(GAOT) to fit diatomic potential energy curves. In order to show this method, we obtain the 
analytical functions of the and Li2 systems using the ab initio energy calculations as well as 
Rydberg trial function. These studies show that the quality of the GAOT fitting is comparable to 
the best optimization techniques recommended to fit diatomic systems. The introduction of this 
new technique is very important because it arises as a new option to fit potential energy surfaces 
for reactive scattering dynamics. 



I. INTRODUCTION 

Genetic algorithms 0, |2i S S B 1^ have been applied successfully in the description of a variety of global mini- 
mization problems as well as have attracted significant attention due to their suitability for large-scale optimization 
problems, specially for those in which a desired global minimum is hidden among many local minima. 

The main object of this paper is to propose a genetic algorithm technique for fitting the potential energy curves 
(PEC) to points obtained by ab-initio calculations. In order to present and to test the method, we reproduce the PEC 
of the following diatomic systems and iz2 using a Rydberg functional form The molecular constants and PEC 
of the and lithium dimer have been of great interest to theoretical chemists, spectroscopists and astrophysicists. The 
i?2^ system is the most simple and probably the most studied molecule. These molecules are relatively small and can 
be treated very accurately. It is not a surprise that they are used for testing and applications of new methodologies 
0. Further, these systems has been intensely investigated in experimental basis. 

Following, we outline how the paper is organized. In Section 2 we present the main characteristics of the GAOT 
used in our calculations. The details of the fitting are shown in Section 3. Our conclusions and detailed comparison 
with other methods are contained in the Section 4. 



II. MODEL 



A. The Problem 



In order to find the PEC of H2'^ and Li2 we want that the GAOT finds a set of parameters [a] — [ai,a2, ■■,am] 
that minimize the mean square deviation of Rydberg functional form. We employed this functional form to adjust 
to ab initio points of these diatomic systems. Given a set of rip ordered pairs (cp, r^) of ab initio points, where is 
the energies and rp is the distance between the nuclei, we want to find a functional form T^([a],r) that minimizes the 
mean square deviation 
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where Cp = V{[a],rp). 

In this work we employ the Rydberg form, y^([a],p), to adjust the potential of diatomic systems through a 
stochastic optimization technique based on GAOT. Our method is general and we can calculate the PES of systems 
with more degrees of freedom. To find the potential form of a given system, we need a set of points that we want to 
adjust. With the total energies as a function of the nuclear configurations we have to choose an appropriate compact 
form to represent them. In this work, the Rydberg form is given by 

m 

V^y{[e.],p) = -De(l + ^a,y)e-»^'' (2) 

The dissociation energy is given by 

De = -V^y{[si],0) (3) 

where p = r — rgg and rgq is the equilibrium bond length of the system. In this work we use m = 3 for the Rydberg 
form. 



B. Codification 



In our GAOT the population is coded in a binary discrete cube named A, with I x m x n bits. The elements of 
A, Uijk, are either or 1, with k integers numbers l<i<n, l<j<m, l<k<n. The label i refers to the 
component i of the gene j of the individual k. Therefore, A represents a population of n individuals, each one of them 
have a genetic code with m genes. Each gene is a binary string with I bits. 

The genetic code of the individual k is given by 

[a]k = [Olfej «2fe, Omfelj (4) 



djk = ^ 2* ^ttijk (5) 

i=l 

is a integer number composed with the binary string aijka2jk--aijk--0'ljk- It is defined on the interval [0,2' — 1]. To 
define the real search space for each parameter, we transform 

- a.fc = ' + (6) 

were a^^ is a real number defined on the interval 

= [a™",a7^"^]. (7) 

Now we define the phenotype of the individual k 

Kp = ^([a]fe, Pkp) = -a5fe(l + aikpkp + a^kPkp^ + a3kPkp^)e~''"""' (8) 

where [a]fe = [aik, a2kj o^akt Ci4k^ o-sklt Pkp = ~ '^•ik and is the interatomic distance of the system. The gene 05;; 
represents the dissociation energy. By this way we have m = m + 2 = 5 independent parameters to optimize. With 
this we define the fitness of a phenotype k 

Fk = Smax — Sk (9) 
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C. Operators 

In our GAOT we use the most common operators: selection, recombination and mutation. The selection operator 
normalize the vector Fj. 

- Tk 

that represents the probability of each individual been selected for a recombination through a roulette spinning. For 
the purpose of this work we selected n/2 individuals (parents) that will generate, through the recombination operator, 
n/2 new individuals (offsprings). So, to make a new generation we joint the n/2 old strings (parents) with an/2 
new strings (offsprings) in order to maintain a population with fixed number n. The recombination operator is a 
cross-over operator that recombine the binary string of each gene j of two random selected individuals to form two 
new individuals. In this work we use a two random point cross-over. 

The mutation operator flip N^nut random selected bits in a population. We choose N^n^t to make the probability 
of change of a given bit equal to 0.01 per cent. So, in a population of I x m x n bits, we make 

(12) 

I X m X n 

where q is the probability of change in one bit. 



D. Linear scaling and elitist strategy 

When we use the GAOT in a minimization procedure we want to find a solution that is a global minimum. In fact, 
this solution can be found by a group of individuals or by all individuals in a population. When all the population 
converges to a single solution, this solution could not be a global minimum, but a local minimum. This is called 
premature convergence and it can be avoided by the linear scaling jj . This procedure enhances the probability that 
several minima will coexists in the population. A general linear scaling in Sk is given by 

S'k = aSk + b. (13) 

In this work we define 

e = ~ (14) 



Sriiax ~ S 
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and if S,nin < e 
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Where Smax and Smin are the worst and the best phenotypes, respectively. S is the mean value of the phenotypes 
of the population. In this way, we maintain the average S'^ = Sk, set S' max = gS if Smax ^ S and S'min = if 
Smin ^ S. 5 is a arbitrary value to control the selective pressure. Expression (14) avoids negative values for 5"^. The 
selection operator normalize the vector F^. — gS — S^. In fact, this linear scaling has lowered the selective pressure 
on the population. Through this, and setting a correct mutation rate, we maintain the variety of the population and 
therefore we avoid the convergence of the population. The elitist strategy consists of copying an arbitrary number 
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FIG. 1: (A) Test function V(x). Frequency of each individual in the population after 250 and 500 generations with (B) linear 
fit on (g=l.f,iVe; = 8) and (C) linear fit off. 

Npi of the best individual on the population in the next generation. It warrants that this individual will not be 
extinguished. Here a example is useful. We want find the minimum of the one dimensional test function 

gg-0.09(x-20)^ _|_ -[^2e~"-09(^-28)^ _|_ 
22g-0.09(x-35)^ _|_ 2g-0-09(a:-50)^ _|_ 
7g-0.2(2;-42.5)" ^ 

at the interval [0, 50]. In Fig. 1 we show the form of this function and the frequency of each individual appearing in 
the population, after 250 and 500 generations. We present the results with S'f. or Sk- It is clear the effect of the linear 
fit in the population dynamics. When we use S*^, after 500 generations, the 880 individuals in the population make a 
signature of the function V(x). 



III. RESULTS 



We run 50 times the algorithm for each case to evaluate the performance of the method. After extensive trials 
of the parameters values we take as parameters for all cases: n=800, Z=30, ra—Q or 8, (7=0.01, A^ei=40 and a 1000 
generations. It should be pointed out that the algorithm is very robust and it works properly with an wide range of 
these parameters. The initial population were always random numbers. 

All the runs produced a subpopulation of acceptable solutions. We defined as acceptable solutions when the root 
mean square (rms) and the difference per point are smaller than 1.0 Kcal/mol. Actually, for each run, we found a great 
number of acceptable solutions. Figure 2 shows the distribution of individuals in the region of interest (acceptable 
solutions) after a 1000 generation of the fifty runs for Li2 PEC. Similar results has been found for H2+. 

Initially we employed the GAOT to obtained the potential curves of ion molecule using a set of highly accurate 
molecular energies (table I of reference 0) and Rydberg function (see equation 2) with r^q and free to change. 
In the Rydberg form r^q and are the genes 04 and 05 respectively. So, we chose their search intervals, 84^ and ^5, 
close to the regions were we have experimental values for these parameters ^ . In Table are shown the parameters 
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500- 




0,0 0,2 0,4 0,6 0,8 1,0 1,2 

Fitness values (Kcal/mol) 



1,4 



1,6 



FIG. 2: Distribution of individuals in the region of interest (acceptable solutions) after a 1000 generation of the fifty runs for 
Li2 PEC. 



obtained by our fitting using the GAOT. The rms deviation value S in this fitting was of the 0.927 x 10^""' Hartree 
(about 0.05810 kcal/mol). The figure 13 shows a comparison between the GAOT and ab initio PEC. From this figure 
one can see that both GAOT and ab initio PEC are in a good agreement. Once a PEC has been fitted to a analytical 
form, both diatomic vibrational energies and spectra can be determined from the radial Schrodinger equation [T(tI |. 
So, to better test our fitting, we resolved the radial Schrodinger equation using our H2'^ PEC through DVR method 
\n\ . The Table lU shows the comparison among the spectra values obtained using the GAOT PEC with the spectra 
of the reference jSj . One can see from these comparison that GAOT spectra are in good agreement with that found 
in the literature. 

To complete our test about of efficiency of the GOAT, we also fit the PEC for the Li2 system from ab initio 
energies. To make this we used the ab initio energies of reference the Rydberg analytical function. We 

obtained good results to this molecule with Veq and De free parameters in the minimization process. In this case, the 
rms deviation value was about 0.3062884602 kcal/mol, Veq = 2.6732799^° and De = 24.44413 kcal/mol. The GAOT 
results obtained for Li2 molecule are in Table IIIII In figure 0| are represented both GAOT and ab initio PEC of the 
Li2 system. From this figure one can see that GAOT and ab initio PEC are in a good agreement. In the Table ITVI 
we compare_the Li2 spectra obtained via GAOT PEC and two others spectra found in the literature, FCIPP [IJ 
and RKR 12], respectively. From this comparison one can note that GAOT spectra are in good agrement with both 
FCIPP and RKR spectra. 



Table I: Parameters obtained of GAOT fitting for in Rydberg form. Energy and nuclear distances are given in atomic 

units. 



Parameters 


H2 + 


ffll 


1.4960937607 




0.5767577799 


aa 


0.21384215222 




1.9970569317 




0.0946586828 
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Table II: H2 spectra given in cm 





GAOT spectra 


Reference ^ 


1 


2190.1104 


2192.022 


2 


4256.9176 


4256.714 


3 


6204.1122 


6198.318 



Table III: Parameters obtained of GAOT fitting for Li2 in a Rydberg form. Energy in Kcal/mol and nuclear distances in 
Angstrons. 



Parameters 


Li2 


ai 


1.91967773 


a2 


1.078125 


as 


0.22248840 




2.69008356 




24.4238281 


S 


0.11981447 




4 6 
Internuclear distance 



FIG. 3: Comparison between the GAOT(solid line) and ab initio{dashed line) PEC of -ff2^ system. 



IV. CONCLUSIONS 



The genetic algorithm is a useful minimization tool. We can use it to obtain the potential curve of diatomic 
systems by adjust ab initio or experimental data. The dynamics of the population in the GAOT given us important 
information about the studied system, in special when we find the distribution of the local minima. Understand this 
dynamics and developing tools and theories to control them is important for employed the GAOT in several classes 
of physical problems. This study was very important, so it enables the GAOT to fit potential energy surface (PES) 
of the bond and scattering systems. 
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Table IV: Li2 spectra given in cm 





LcAU i spectra 


r UIFFI i/] 


Kls.Kll/| 


1 


347.42 


346.05 


346.46 


2 


688.69 


686.65 


687.86 


3 


1023.79 


1021.71 


1024.08 


4 


1352.65 


1351.15 


1355.01 


5 


1675.21 


1674.88 


1680.54 


6 


1991.43 


1992.81 


2000.56 


7 


2301.24 


2304.85 


2314.95 


8 


2604.59 


2610.92 


2623.58 


9 


2901.42 


2910.90 


2926.35 


10 


3191.66 


3204.70 


3223.11 




2 4 6 8 

Internuclear distance 



FIG. 4: Comparison between the GAOT(solid line) and ah mrfio(daslied line) PEC of Li2 system. 
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